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In almost all of the studies devoted to the time delay and the frequency shift of light, the calcu- 
. . . lations are based on the integration of the null geodesic equations. However, the above-mentioned 

C^^ ' effects can be calculated without integrating the geodesic equations if one is able to determine the 

^^ , bifunction Q,{xA,oi:B) giving half the squared geodesic distance between two points xa and xb (this 

^^ . bifunction may be called Synge's world function). In this lecture, Q{xa, xb) is determined up to the 

Cn ' order 1/c? within the framework of the PPN formalism. The case of a stationary gravitational field 

generated by an isolated, slowly rotating axisymmetric body is studied in detail. The calculation 
^^ ■ of the time delay and the frequency shift is carried out up to the order 1/c^. Explicit formulae are 

Q_^ ' obtained for the contributions of the mass, of the quadrupole moment and of the internal angular 

momentum when the only post-Newtonian parameters different from zero are jS and 7. It is shown 
that the frequency shift induced by the mass quadrupole moment of the Earth at the order 1/c? will 
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amount to 10 in spatial experiments like the ESA's Atomic Clock Ensemble in Space mission. 
Other contributions are briefly discussed. 



I. INTRODUCTION 



^ . A lot of fundamental tests of gravitational theories rest on highly precise measurements of the travel time and/or 

T^ ' the frequency shift of electromagnetic signals propagating through the gravitational field of the Solar System. In 
C^ . practically all of the previous studies, the explicit expressions of such travel times and frequency shifts as predicted 
^D ' by various metric theories of gravity are derived from an integration of the null geodesic differential equations. This 
^—^ \ method works quite well within the first post-Minkowskian approximation, as it is shown by the results obtained, 
e.g., in [1-5]. Of course, it works also within the post-Newtonian approximation, especially in the case of a static, 
spherically symmetric space-time treated up to order l/c' [6, 7]. However, the solution of the geodesic equations 
l^ . requires heavy calculations when one has to take into account the presence of mass multipoles in the field or the tidal 
effects due to the planetary motions, and the calculations become quite complicated in the post-post-Minkowskian 
approximation [8] , especially in the dynamical case [9] . 

The aim of this lecture is to present a quite different procedure recently developed by two of us. Based on Synge's 
world function [10], this procedure avoids the integration of the null geodesic equations and is particularly convenient 
J^ ' for determining the light rays which connect an emitter and a receiver having specified spatial locations at a finite 
distance. Thus, we are able to extend the previous calculations of the time delay and of the frequency shift up to the 
order l/c"*. As a consequence, it is now possible to predict the time/frequency transfers in the vicinity of the Earth at 
a level of accuracy which amounts to 10~^^ in fractional frequency. This level of accuracy is expected to be reached 
in the foreseeable future with optical atomic clocks [11]. 

The plan of the lecture is as follows. First, in Sect. II, the definition of the time transfer functions are given and the 
invariant expression of the frequency shift is recalled. It is shown that explicit expressions of the frequency shift can 
be derived when the time transfer functions are known. In Sect. Ill, the relevant properties of Synge's world function 
are recalled. In Sect. IV, the general expressions of the world function and of the time transfer functions are obtained 
within the Nordtvedt-Will parametrized post-Newtonian (PPN) formalism. In Sect. V, the case of a stationary field 
generated by an isolated, slowly rotating axisymmetric body is analyzed in detail. It is shown that the contributions 
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of the mass and spin multipoles can be obtained by straightforward derivations of a single function. Retaining only 
the terms due to the mass M, to the quadrupolc moment J2 and to the intrinsic angular momentum S of the rotating 
body, explicit expansions of the world function and of the time transfer function are derived up to the order l/c" and 
1/c*, respectively. The same formalism yields the vectors tangent to the light ray at the emitter and at the receiver 
up to the order \/c^ . In Sect. VI, the frequency shift is developed up to the order 1/c^ on the assumption that (3 and 
7 are the only non vanishing post-Newtonian parameters. Explicit expressions are obtained for the contributions of 
J2 and S. Numerical estimates are given for the ESA's Atomic Clock Ensemble in Space (ACES) mission [12, 13]. 
Concluding remarks are given in Sect. VII. 

Equivalent results formulated with slightly different notations may be found in [14] and an extension of the method 
to the general post-Minkowskian approximation is given in [15]. 

Notations. - In this work, G is the Newtonian gravitational constant and c is the speed of light in a vacuum. The 

Lorentzian metric of space-time is denoted by g. The signature adopted for g is (H ). We suppose that the 

space-time is covered by one global coordinate system (x^) = {x^^x), where x" — ct, t being a time coordinate, and 
X = (x*), the X* being quasi Cartesian coordinates. We choose coordinates x* so that the curves of equations x* — 
const are timclike. This choice means that goo > everywhere. We employ the vector notation a in order to denote 
either {a^, a^, a^} — {a*} or {ai, a2, 03} — {ui}. Considering two such quantities a and b with for instance a = {a'}, 
we use a • b to denote a*6* if b = {&*} or a^bi ii b — {bi} (the Einstein convention on the repeated indices is used). 
The quantity | a \ stands for the ordinary Euclidean norm of a. 

II. TIME TRANSFER FUNCTIONS, TIME DELAY AND FREQUENCY SHIFT 

We consider here electromagnetic signals propagating through a vacuum between an emitter A and a receiver B. We 
suppose that these signals may be assimilated to light rays travelling along null geodesies of the metric (geometic optics 
approximation). We call xa the point of emission by A and Xb the point of reception by B. We put xa = {ctA, xa) 
and Xb = {cIbiXb)- We assume that there do not exist two distinct null geodesies starting from xa and intersecting 
the world line of B. These assumptions are clearly satisfied in all experiments currently envisaged in the Solar System. 

1. Time transfer functions and time delay- The quantity t^ ^ ^A is the (coordinate) travel time of the signal. 
Upon the above mentioned assumptions, tB — tA may be considered either as a function of the instant of emission tA 
and of xai xb, or as a function of the instant of reception ts and of xa and xb- So, we can in general define two 
distinct (coordinate) time transfer functions, 71 and Tr by putting : 

tB -tA = %{tA, Xa, Xb) = Tr{tB,XA, Xb) ■ (1) 

We call Te the emission time transfer function and 7^ the reception time transfer function. As we shall see below, 
the main problem will consist in determining explicitly these functions when the metric is given. Of course, it is in 
principle sufficient to determine one of these functions. 
We shall put 

Rab =\xb-xa\ ■ (2) 

throughout this work. The time delay is then defined as tB — tA — Rab/c. It is well known that this quantity is > 
in Schwarzschild space-time, which explains its designation [16]. 

2. Frequency shift- Denote by u^ and u^ the unit 4- velocity vectors of the emitter at xa and of the receiver 
at Xb, respectively. Let Tab be the null geodesic path connecting xa and xb, described by parametric equations 
x" = x"(C), C being an affine parameter. Denote by l'^ the vector tangent to Tab defined as 

l'-^- (3) 

dC 

Let va be the frequency of the signal emitted at xa as measured by a clock comoving with A and ub be the frequency 
of the same signal received at Xb as measured by a clock comoving with B. The ratio va/vb is given by the well-known 
formula [10] 

VA_ _ u'X{1^)a , . 

l^B Wg{l^)B 

Since it is assumed that the emission and reception points are connected by a single null geodesic, it is clear that 
{1^)a and (?^)b may be considered cither as functions of the instant of emission tA and of xa, xb, or as functions of 
the instant of reception Ib and of xa and xb- Therefore, we may write 

^ Ne{uA,UB;tA,XA,XB) ^ Afr{uA,UB;tB, Xa,Xb) ■ (5) 

l^B 



Denote by va = {dx/dt)A and vb = {dx/dt)B the coordinate velocities of the observers at xa and XB^ respectively: 

(6) 



VA = 



dx 



It is easy to see that the formula (4) may be written as 

i^A _ Ma ('o)a Qa 



Vb = 



1, 



dx 



VB u% (?o)b qB 
where I a and Ib are the quantities defined as 
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1 + -Ib ■ Vb , 
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(7) 



(8) 



It is immediately deduced from Eq. (7) that an explicit expression of A/'e (resp. A/'r) can be derived when the time 
transfer function 7^ (resp. %•) is known. Indeed, one has the following theorem [15]. 

Theorem. Consider a signal emitted at point xa — {ctA,XA) and received at point xb = {ctB,XB)- Denote by /^ 
the vector dx^ / dQ tangent to the null geodesic at point x{Q), Q being any affine parameter, and put 



h 



Then, one has relations as follow at xa and at xb 
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(9) 

(10) 

(11) 
(12) 



where %. andTr are taken at (i^i, a;^, aj^) and {tB,XA,XB), respectively. 

This theorem may be straightforwardly deduced from a fundamental property of the world function that we intro- 
duce in the following section. 

Case of a stationary space-time. In a stationary space-time, we can choose coordinates (x^) such that the metric 
does not depend on x'^ . Then, the travel time of the signal only depends on xa, xb. This means that Eq. (1) reduces 
to a single relation of the form 



tB -tA = T{xa,xb) ■ 
It immediately follows from Eqs. (10) and (11) that 
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-c^-r'^ixA,XB) 

OXg 
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As a consequence, the formula (7) reduces now to 



VA 

Vb 



i\ l + VA-V^^T 

1% 1-Vb- Va^flT 



(13) 

(14) 
(15) 
(16) 

(17) 



where V^/ denotes the usual gradient operator acting on f{x). 

It is worthy of note that (1, {(^i)^}) and (1, {{li)B}) constitute a set of covariant components of the vector tangent 
to the light ray at xa and Xb, respectively. This tangent vector corresponds to the afhne parameter chosen so that 
{Io)a = (/o)b = 1. 



III. THE WORLD FUNCTION AND ITS POST-NEWTONIAN LIMIT 
A. Definition and fundamental properties 

For a moment, consider xa and Xb as arbitrary points. We assume that there exists one and only one geodesic 
path, say Tab, which links these two points. This assumption means that point xb belongs to the normal convex 
neighbourhood [17] of point xa (and conversely that xa belongs to the normal convex neighbourhood of point xb)- 
The world function is the two-point function fl(xA,XB) defined by 

n{xA,XB) ^ -eAsisAB]'^ , (18) 

where sab is the geodesic distance between xa and xb, namely 



SAB = / ^Jg^,udx^'dx'' (19) 

JTab 

and eAB ~ l^O,— 1 according as Tab is a timelike, a null or a spacelike geodesic. An elementary calculation shows 
that U,{xA, Xb) may be written in any case as [10] 



1 [^ dx^ dx^ 

n{xA,XB)^ -J g^^{x"(X)) — —dX, 



(20) 



where the integral is taken along Tab , ^ denoting the unique affine parameter along Tab which fulfills the boundary 
conditions A^ = and Xb — 1- 

It follows from Eqs. (18) or (20) that the world function fl(xA,XB) is unchanged if we perform any admissible 
coordinate transformation. 

The utility of the world function for our purpose comes from the following properties [10, 15]. 

i) The vectors (dx" /dX)A and {dx°'/dX)B tangent to the geodesic Tab respectively at xa and xb are given by 

As a consequence, if f2(x^, xb) is explicitly known, the determination of these vectors does not require the integration 
of the differential equations of the geodesic. 

a) Two points XA and xb are linked by a null geodesic if and only if the condition 

n{xA,XB)^0 (22) 

is fulfilled. Thus, Vt{xA, x) — Q \& the equation of the null cone C{xa) at xa- 

Consequently, if the bifunction U,{xa,xb) is explicitly known, it is in principle possible to determine the emission 
time transfer function 7^ by solving the equation 

VL{ctA,XA,ctB,XB) = Q (23) 

for Ib- It must be pointed out, however, that solving Eq. (23) for Ib yields two distinct solutions t^ and t^ since 
the timelike curve x^ = x^ cuts the hght cone C{xa) at two points x^ and a;^, x^ being in the future of x^. Since 
we regard xa as the point of emission of the signal and Xb as the point of reception, we shall exclusively focus our 
attention on the determination of ig — Ia (clearly, the determination of t^ — tA comes within the same methodology). 
For the sake of brevity, we shall henceforth write ts instead of t~^ . 

Of course, solving Eq. (23) for tA yields the reception time transfer function %. 

Generally, extracting the time transfer fonctions from Eq. (23), next using Eqs. (10) or (11) will be more straight- 
forward than deriving the vectors tangent at xa and xb from (21), next imposing the constraint (22). 

To finish, note that the theorem stated in Sect. II is easily deduced from the identities 

n{ctA,XA,ctA + C%{tA,XA,XB),XB) =0 (24) 

and 

ft{ctB - cTr{tB,XA,XB),XA,ctB,XB) = 0- (25) 



B. General expression of the world function in the post- Newtonian limit 

We assume that the metric may be written as 



9^„y = V^ii' + hf_,^ (26) 

7]^,^ =diag(l,-l,- 
equations x" — x?qJX), with 



throughout space-time, with rjf^i, =diag(l, — 1, — 1, — 1). Let Fj^^ be the straight hne defined by the parametric 



x1o^{\)^{x%-x'X)\ + x'X, 0<A<1. (27) 

With this definition, the parametric equations of the geodesic Tab connecting xa and xb may be written in the form 

x"{\)=xf,^{\)+X"{X), 0<A<1, (28) 

where the quantities X"{X) satisfy the boundary conditions 

X"(0)=0, X"(1) = 0. (29) 

Inserting Eq. (26) and dx'^{X)/dX = Xg — Xj^ + dX'^{\)/d\ in Eq. (20), then developing and noting that 

1 dX'^ 

rj^,{x%-x'X)—-dX^O (30) 

"'■^ 



by virtue of Eq. (29), we find the rigorous formula 
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nixA,XB) - n'^'''>{xA,XB) + -{x''s^A)i=^B-^A) K.{x'^{X))dX 

^ Jo 



-I 

'2 Jo 



5..(^"(A))^^ + 2(x^ - x^) V(^"(A))'^^' 



dX, (31) 



dX dX ' " -^,-^^y- vv. ^^ 

where the integrals are taken over Tab and Vv-^'{xA^ xb) is the world function in Minkowski space-time 

f](")(x^,XB) - \vAx^B - <Wb - ^a) ■ (32) 

Henceforth, we shall only consider weak gravitational fields generated by self-gravitating extended bodies within 
the slow-motion, post-Newtonian approximation. So, we assume that the potentials ft,^^ may be expanded as follows 

hoo = ^hl^^ + ^/i^o' + 0(6) , 

1 -,(3) 



ho^ = ^K^ + 0(5) , (33) 

From these expansions and from the Euler-Lagrange equations satisfied by any geodesic curve, namely 

d { dx/^X 1^ , dxf' dx" 

Tx['-'^)^2^-^'^^^' (^'^ 

it results that X^( A) = 0(2) and that dxi^/dX = a;^ -.t^ + 0(2). As a consequence, /i^^(a:;"(A)) = h^^{x?^AX)) + 0{'i) 
and the third and fourth terms in the r.h.s. of Eq. (31) arc of order 1/c^. These features result in an expression for 
^{xa,xb) as follows 

VL{xA.XB)^n^''\xA,XB) + n^''''HxA.XB) + 0{A), (35) 

where 

n''P''\xA,XB) = ^(4 -0'^'4o^(4o)(A))dA 

+ ^{^B - ^a){^'b A) t hif(x1,^{X))dX 

+ ^{x% x\){x^B A) I 4?(x^o)(A))rfA, (36) 



c 







the integral being now taken over the Hne FJ^^ defined by Eq. (27). 

The formulae (35) and (36) yield the general expression of the world function up to the order 1/c^ within the 
framework of the 1 PN approximation. We shall see in the next subsection that this approximation is sufficient to 
determine the time transfer functions up to the order 1/c^. It is worthy of note that the method used above would 
as well lead to the expression of the world function in the linearized weak-field limit previously found by Synge [10]. 

C. Time transfer functions at the order 1/c* 

Suppose that xb is the point of reception of a signal emitted at xa- Taking Eq. (35) into account, Eq. (22) may 
be written in the form 

n^'>\xA,XB) + n^''''HxA,XB) = 0(4) , 

which implies the relation 

tB-tA^ '-Rab - "^""nc^.^a..,ct.,..) ^ ^^^ ^3^^ 

c cRab 

Using iteratively this relation, we find for the emission time transfer function 

n-U \ Id ^^^^\ctA,XA,ctA + RAB,XB) , „,.. .„^. 

%(tA,XA,XB) = -Rab 7^ hO(5). (38) 

c cKab 

and for the reception time transfer function 

n-U \ ^D n'^^^^ctB-RAB,XA,ctB,XB) ,^..^ .on^ 

%.{tB,XA,XB) = -Rab 5 hO(5). (39) 

c cRab 

These last formulae show that the time transfer functions can be explicitly calculated up to the order l/c"* when 
Q,\P^)(^xa,xb) is known. This fundamental result will be exploited in the following sections. 

It is worthy of note that a comparison of Eqs. (38) and (39) immediately gives the following relations : 

%{tB,XA,XB) ^ rJtB - ^,XA,XBJ + 0(5) (40) 

and conversely 

%{tA,XA,XB) =Tr(tA + -^^,Xa,XbJ +0(5). (41) 



The quantity H,^^^' {ctA, xa, cIa + Rab, xb) in (38) may be written in an integral form by using Eq. (36), in which 

(0)( 



Rab and Rab^ + c^a are substituted for x^ — x\ and for x''LA\), respectively. As a consequence 



%{tA, XA, Xb) = -Rab \^- ^ 



h^^JizUX)) + hf{zUX))N^N^ + -4?«(A))iV^ 
■' c 



d\ } + 0(5) , (42) 



the integral being taken over curve T^^ defined by the parametric equations x" = z"(A), where 

zl{\) = RAB\ + ctA, z\{\)=RabN'\ + x\, 0<A<1, (43) 

with 

Rab^Rab\, N^^'^M^^, (44) 

Rab 

We note that T a^ is a null geodesic path of Minkowski metric from xa, having the above-defined quantities TV* as 
direction cosines. 

A similar reasoning leads to an expression as follows for 7^ 



%.{tB,XA,XB) = -Rab ] ^ ^ T^ / 



hg\z^{X)) + hf(z'l{\))Nm= + -4f (zf (A))7V* 
■^ c 



d\ \ + 0(5) , (45) 



the integral being now taken over curve Ta^ defined by the parametric equations x" = z"(A), where 

z°(A) = -i^^sA + cis, zI{\) = -RabN^\ + x^b, 0<A<1. (46) 

Curve T AB ^^ ^ ^^^^ geodesic path of Minkowski metric arriving at xb and having A^* as direction cosines. 



IV. WORLD FUNCTION AND TIME TRANSFER FUNCTIONS WITHIN THE NORDTVEDT-WILL 

PPN FORMALISM 

A. Metric in the 1 PN approximation 



In this section, we use the Nordvedt-WiU post-Newtonian forniahsni involving ten parameters /3, 7, ^, ai, . . ., ^4 
[18]. We introduce shghtly modified notations in order to be closed of the formalism recently proposed by Klioner 
and Soffel [20] as an extension of the post-Newtonian framework elaborated by Damour, Soffel and Xu [21] for general 
relativity. In particular, we denote by Vr the velocity of the center of mass O relative to the universe rest frame [27]. 

Although our method is not confined to any particular assumption on the matter, we suppose here that each source 
of the field is described by the energy-momentum tensor of a perfect fluid 



rpflV 



pc 



i + l(n + ^ 
c V pj \ 



u^u" - pg"" , 



where p is the rest mass density, 11 is the specific energy density (ratio of internal energy density to rest mass density), 
p is the pressure and u^ is the unit 4-velocity of the fluid. In this section and in the following one, v is the coordinate 
velocity dx/dt of an element of the fluid. We introduce the conserved mass density p* given by 

^ y + 3"/U] +0(4) , 



P =PV-9u 



1 



where g = dct{g^i,) and U is the Newtonian- like potential 



U{x",x)=G 



P*{x\x') 
\x — x'\ 



dV. 



(47) 



(48) 



In order to obtain a more simple form than the usual one for the potentials h^i, we suppose that the chosen (x^) 
are related to a standard post- Newtonian gauge (x^) by the transformation 



1 



a;0 - xO + ^ [(1 + 2i + a2- Ci)dtX - 2a2'yr • Vx] , 
where x is the superpotential defined by 

xix",x) = ^Gjp*{x^,x')\x-x'\d^x'. 
Moreover, we define p by 



(49) 
(50) 



P = P 



1 + ^(27 + 1 - 2C + ^3 + Ci)^ + (1 - 2/3 + e + C2)^ + (1 + C3)^ + (37 - 2^ + Ki)^ 



p^c- 



-(ai - aa)-! - -[ax - 203)^-^— 



0(4) 



Then, the post-Newtonian potentials read 

2/3 



^oo 



2e 



2a2 



^ + ^w' + ^<\>w + ^(Ci - 200. - -^v\vld.,x + 0(6) 



h = {hoi\ = 



1 \ 1 

7 + 1 + -ai 1 w + -aiwVr 



0(5), 



where 



27 

Kj = — 5"^% + 0(4) 



w{x\x)^G I ^^^^d^x' 



1 



(j)w{x^,x) = G 



2 ^(1 + 2^ + a2 - Ci)dttX - "^oi-iVr ■ V{dtx)\ , 
2 ( p*{x\x')p*{x\x"){x-x') 



X — X' \ c^ 



\X — X 



/|3 



X — X" 



d^x'd^x". 



Kix^x) . G I p-(^°.^ob(^°.^;H^-^o]-,3,, 



XX 

^1 



'13 






(51) 

(52) 
(53) 
(54) 

(55) 
(56) 
(57) 
(58) 



B. Determination of the world function and of the time transfer functions 

For the post-Newtonian metric given by Eqs. (52)-(58), it follows from Eq. (36) that fl{xA,XB) may be written up 
to the order 1/c^ in the form given by Eq. (35) with 



n^''''\xA,XB) = nl''''HxA,XB) + n(P''\xA,XB) + n[!:''HxA,XB). 



where 



f^(f^)(xA,XB) = --[(4-4)'+7i?iB] / «^(xfo)(A))ciA, 
c Jo 

r!(f ^) {xA,XB)^^(l+l + \a^ {x% - x\)Rab ■ j w{xf,^ {X))dX , 

1 r^ 

nif/^XA^XB) = Tr^ai{x% - x"^){Rab ■ Vr) / «;(xfo)(A))dA , 
^c Jq 



the integrals being calculated along the curve defined by Eq. (27). 

The emission time transfer function is easily obtained by using Eqs. (38) or (42). We get 



11 f^ 

Te{tA,XA,XB) = -RaB + —{i + 1)RaB I 
C C Jq 



tRab ■ 



w{z^{X))dX 



(7 + l + -ai)/ w{zl{X))dX + -aiVr ii;(z^(A))dA 



the integral being evaluated along the curve fJ^^ defined by Eq. (43). 
The reception time transfer function is given by 



1 



1 



TritB,XA,XB) = -RaB + ^il + 1)RaB w(z^(A))dA 



c 
2 



tRab 



(7 + 1 + iai) / w{z'^{X))dX + -aiVr j w{z'l{\))dX 



0(5), 



0(5)^ 



(59) 

(60) 
(61) 
(62) 



(63) 



(64) 



the integral being evaluated along the curve F)^^ defined by Eq. (46). 

Let us emphasize that, since w = U + 0(2), w may be replaced by the Newtonian- like potential U in Eqs. (60)-(63). 



C. Case of an stationary source 

In what follows, we suppose that the gravitational field is generated by a single stationary source. Then, dtX — 
and the potentials w and w do not depend on time. In this case, the integration involved in Eqs. (60)-(62) can be 
performed by a method due to Buchdahl [19]. Introducing the auxiliary variables y^ = xa — x' and ys = Xb — x' , 
and replacing in Eq. (27) the parameter A by u = A — 1/2, a straightforward calculation yields 



w(a;(o)(A))(iA — G p{x')F{x' ,XA,XB)d^x' , 



w{x(^Q-^{Xj)dX = G p*{x')v{x')F{x',XA,XB)d^x' 



where the kernel function F(x' , xa, xb) has the expression 

.1/2 



F{x',xa,xb) = / 

J-1/ 



du 



1/2 I (ye - yA)u + ^(ys + va) 

Noting that ys — yA = Rab, which implies that | ys — yA \= Rab, we find 



Rai 



\x-xa\ + \x-xb\ -Rai 



(65) 
(66) 

(67) 

(68) 



Inserting Eqs. (65), (66) and (68) in Eqs. (60)-(62) and in Eq. (63) will enable one to obtain quite elegant 
expressions for D,^-^^\xa,xb) and for T{xa,xb), respectively. 



V. ISOLATED, SLOWLY ROTATING AXISYMMETRIC BODY 

Henceforth, we suppose that the hght is propagating in the gravitational field of an isolated, slowly rotating 
axisymmetric body. The gravitational field is assumed to be stationary. The main purpose of this section is to 
determine the infiuence of the mass and spin multipole moments of the rotating body on the coordinate time transfer 
and on the direction of light rays. From these results, it will be possible to obtain a relativistic modelling of the 
one-way time transfers and frequency shifts up to the order 1/c^ in a geocentric non rotating frame. 

Since we treat the case of a body located very far from the other bodies of the universe, the global coordinate 
system (a;^) used until now can be considered as a local (i.e. geocentric) one. So, in agreement with the UAI/UGG 
Resolution Bl (2000) [22], we shall henceforth denote by W and W the quantities w and w respectively defined by 
Eqs. (55) and (58) and we shall denote by Gfj,i, the components of the metric. However, we shall continue here with 
using lower case letters for the geocentric coordinates in order to avoid too heavy notations. 

The center of mass O of the rotating body being taken as the origin of the quasi Cartesian coordinates (a;), we 
choose the axis of symmetry as the x'^-axis. We assume that the body is rotating about Oa;"^ with a constant angular 
velocity a;, so that 



v(x) 



(69) 



In what follows, we put r =| a; |, r^ —\ xa \ and rg =| xa \- We call 6 the angle between x and Oa;"^. We consider 
only the case where all points of the segment joining x^ and Xb are outside the body. We denote by r^ the radius 
of the smallest sphere centered on O and containing the body (for celestial bodies, r^ is the equatorial radius). In 
this section, we assume the convergence of the multipole expansions formally derived below at any point outside the 
body, even if r < re . 



A. Multipole expansions of W and W 



According to Eqs. (55), (58) and (69), the gravitational potentials W and W obey the equations 

V^iy ^ -inGp, V^W = -AttGp*uj x x . 



(70) 



It follows from Eq. (70) that the potential M^ is a harmonic function outside the rotating body. As a consequence, 
W may be expanded in a multipole series of the form 



W{x) 



CM 



l-£j„ 



P„ (cos i 



In this expansion, P„ is the Legendre polynomial of degree n and the quantities M, J2, 
generalized Blanchet-Damour mass multipole moments in general relativity [23]. 
For the sake of simplicity, put 



(71) 
J„, . . . correspond to the 



z — X 



Taking into account the identity 



9" /I 
9z" \ r 



(-l)"n! 



^l-\-n 



Pn{z/r), z^x^ 



r / ^ 



it may be seen that 

W{x) = GM 

Substituting for W from Eq. (72) into Eq. (70) yields an expansion for p as follows 

p{x) = M 



r ^ — ' n! ^ dz^' \ r 

n—2 






(72) 



(73) 



where 6^^'{x) is the Dirac distribution supported by the origin O. This expansion of p in a multipole series will be 
exploited in the next subsection. 



Now, substituting Eq. (69) into Eq. (58) yields for the vector potential W 

W(x) = G I ■ 

J \x-x' 

It is possible to show that this vector may be written as 



p*{ x')uixx' ^^^, ^ 



W = --ujx VV, 



10 



(74) 



(75) 



where V is an axisymmetric function satisfying the Laplace equation V^V = outside the body. Consequently, we 
can expand V in a series of the form 



V{x) = ^ 
r 



X — ^ / r \ """ 



n=l 



(76) 



where / and each Kn are constants. Substituting for V from Eq. (76) into Eq. (75), and then using the identity 

{n + l)P„(z/r) + {z/r)P'^{z/r) = PUii^/r) , 
we find an expansion for W as follows 

GIus X X 

1 - 



W{x) 



2r3 



1-Y.Kn 



' PUiicosO) 



(77) 



which coincides with a result previously obtained by one of us [24] . This coincidence shows that I is the moment of 
inertia of the body about the z-axis. Thus, the quantity S ~ lui is the intrinsic angular momentum of the rotating 
body. The coefhcients Kn are completely determined by the density distribution p* and by the shape of the body 
[24, 25]. The quantities /, Ki^ K2, ■■■Kn, ■■■ correspond to the Blanchet-Damour spin multipoles in the special case of 
a stationary axisymmetric gravitational field. 
Equation (77) may also be written as 



W{x) 



1 



GSxV 



-. 00 



i-ir,^. a- /I 



r •' — ' n! '^ 9z" \ r 

n— 1 



Consequently, the density of mass current can be expanded in the multipole series 



p*{x){u3 X x) 



--5 XV 



(-1) 



Qr. 



^^^^(-)-E^^"<'^^'^H- 



(78) 



(79) 



This expansion may be compared with the expansion of p given by Eq. (73). 



B. Multipole structure of the w^orld function 

The function fl^^^\xA,XB) is determined by Eqs. (59)-(62) where w and w are respectively replaced by W and 
W. The integrals involved in the r.h.s. of Eqs. (59)-(62) are given by Eqs. (65) and (66). Substituting Eq. (73) into 
Eq. (65) and using the properties of the Dirac distribution, we obtain 



W{x^o)W)dX = GM 






ri=2 



F{x,xa,xb] 



x=0 



Similarly, substituting Eq. (79) into Eq. (66), we get [28] 



W{x^o)W)dX^-GSxV 






F{x,xa,xb) 



(80) 



(81) 



These formulae show that the multipole expansion of il,^-^^^ (xa , xb) can be thoroughly calculated by straightforward 
differentiations of the kernel function F{x,xatXb) given by Eq. (68). They constitute an essential result, since they 
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give an algorithmic procedure for determining the niultipolc expansions of the time transfer function and of the 
frequency shift in a stationary axisymmetric field (see also Rcf. [2]). 

In order to obtain explicit formulae, we shall only retain the contributions due to M, J2 and S in the expansion 



yielding ri\y and fi 



JPN) 



Then, denoting the unit vector along the z-axis by k and noting that S = Sk, we get for 



JPN) 



^w (xa^xb) 



GM {x% - x^)^ + jR\g ^^ rrA + rB±RAB\ 
Rab KrA + rs- Rab ) 



+ 



2GM 
GM 



-J2r: 



2 i^B 



„0 \2 



J2r: 



[{rA + TB 

(4 



'''''- ir A +rB) 



|2 - /?2 
' ^AB 



rx, ' X j\ fv ' X J^ 



-/R 



AB 



C^ 



{rA+rB?-R\B 



{k X x^ 



rA rB 

(k X xbY 



and for ilL \xa,xb) 



'VK 



Q. 



(PN) 
W 



{xa,xb) 



7 + 1 + -ai 



2GS 



{x"b 



rA+rB fc ■ (XA X xi 
rArB (rA+rB) 



/?2 
^AB 



+ 



(82) 



(83) 



Finally, owing to the limit | ai |< 4 x 10 ^ furnished in [18], we shall henceforth neglect all the multipole contributions 
in Qi,^ '{xa,xb)- Thus, we get 



GM 



niPN){xA.XB)=a^^-^{x%-x\) 



X Rab -Vr , f rA + rB + Rai 



■In 



2(,3 ^--t! -Aj j^^^ — 1 ^ + rB - Rab 



(84) 



In this section and in the following one, the symbol + • • • stands for the contributions of higher multipole moments 
which are neglected. For the sake of brevity, when + • • • is used, we systematically omit to mention the symbol 0{n) 
which stands for the neglected post-Newtonian terms. 



C. Time transfer function up to the order 1/c* 



In what follows, we put 



riA = 



XA 

rA ' 



riB = 



Xb_ 

rB 



and 



Nab = {N'j 



Xb -xa 
Rab 



Furthermore, we use systematically the identity 

[rA + rsf - Rab = 2rAri3(l + ua ■ ub) 



(85) 



(86) 



(87) 



By substituting Rab for x^ — x^ into Eqs. (82)- (84) and inserting the corresponding expression of H,^^^' into Eq. 
(38), we get an expression for the time transfer function as follows 



where 



T{xa,xb) = -Rab +Tm{xa,xb) +Tj^{xa,xb) +Ts{xa,xb) +%^{xa,xb) 
c 



rr ( X . ^lA^^l ('<-A + rB+RAB 

Tm{xa,xb) = (7 + 1)^^ In 



rA + rB - Rai 



rr , N l+\GM ^ rl Rab 

Tj2[Xa,Xb) = -:^J2- 



2 c3 rArB l + UA-TlB 



1 1 \ (fc • Ua + k ■ Ub) 



rA rB J l + UA-TlB 

1 - (fc • UAf 1 - (fc • nB)2 



rA 



rB 



(88) 



(89) 



(90) 



%^{xa,xb) = 



-(7 + 1 
GM 



-ai 



GS 



ai—-^{NAB ■Vr)ln 



1 1 \ fc • {ua X Ub) 

TA TB ) l + UA-riB 

TA + rB + Rab ' 



rA + rB - Rai 
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(91) 
(92) 



The time transfer is thus explicitly determined up to the order 1/c^. The term of order 1/c^ given by Eq. (89) is 
the well-known Shapiro time delay [16]. Equations (90) and (91) extend results previously found for "f — 1 and ai = 
[1]. However, our derivation is more straightforward and yields formulae which are more convenient to calculate the 
frequency shifts. As a final remark, it is worthy of note that Tm and Tj^ are symmetric in {xa,Xb), while Tg and 
7i^ are antisymmetric in {xa,xb)- 



D. Directions of light rays at xa and xb up to the order l/c^ 

In order to determine the vectors tangent to the ray path at xa and xb, we use Eqs. (14) and (15) where T is 
replaced by the expression given by Eqs. (88)-(92). It is clear that Ia and Ib may be written as 



Ia = -Nab + K(xa-,xb); 
Ib = -Nab + K{xa,xb)- 



(93) 
(94) 



where Ae and A^ are perturbation terms due to Tm, Tj^, Ts, Tk„, ■■■ For the expansion of T given by Eqs. (88)-(92), 
we find 



Xe{xA,XB) = -Xm{xb,xa) - Xj^ixB , Xa) + Xs{xb , xa) +XvAxb,xa) 



(95) 



XrixA,XB) = XMixA,XB) +X,j^{xa,Xb) + Xs{xa,Xb) +X.u^{xa,Xb) H , (96) 

where Am, Xj^, Xs and A^,^ stand for the contributions of TJvf, Tj.-^, Ts and 7^^, respectively. We get from Eq. (89) 

1 /-- Rab 



..GM ( 1 1 
Xm{xa,xb) = -(7 + 1)^5- [z— + ZT' 



TA TB I l + riA-nB 



Nai 



TA + rB 



-riB 



(97) 



From Eq. (90), we get 



GM f 1 



Xj^{xa,xb) = (7 + 1)^^ \ J; 



1 



1 



X {Nab 



TA tbJ tatb {l + UA-riB) 
(k-UA + k ■ub)'^ fvA rs 1 



TB rA 2 2 



l + UA-riB 

1 rArB fl-{k-nAY , I- {k-UB^X (ta , tb 



2rA+rB 
Rab 



rA 



TB 



1 h 1 - riA • "B 

TB TA 



-riB- 



TA + rB 



{k-UA + k-UBY (vA rB 3 1 

\ \ ua ■ riB 

l + UA-TiB \rB rA 2 2 



1 ri -iii, „ ^2l rAi2 + nA-nB)+rB 
- [1 - 3(fc • ub) \ — 

1 /l-(fc-nA)' 2{k ■ nA){k ■ tib) 
'oirA + tb) 

2 \ TA TB 



+k 



Rai 



TB 



(fc -Ua) + [k ■ Ub) 



rA{2 + nA -Ub) +rB 



TA + rB 



(98) 
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From Eqs. (91) and (92), we derive the other contributions that are not neglected here : 



>^sixA,XB) = ( 7 + 1 + -a 



GS 



1 



1 



c^tb \rA tb I l + rbA-ns 



1 



X -^ fe X Ua — 



k ■ (UA X Ub) 
l + UA-riB 



rAi2 + nA -tib) +rB 

UA ^ ; riB 



TA + rB 



>^VrixA,XB) = ai 



GM 
~2^ 



Vr - {vr ■ Nab) Nab ^^ [ta + rB + Rab 



Rab 

{vr ■ Nab) 



(99) 



J A + rB - Rab, 

l + nA-riB \rA ^b J \ ta + tb 

We note that the mass and the quadrupole moment yield contributions of order 1/c^, while the intrinsic angular 



-Ub 



(100) 



momentum and the velocity relative to the universe rest frame yield contributions of order 1/c . 



E. Sagnac terms in the time transfer function 



In experiments like ACES Mission, recording the time of emission Ia will be more practical than recording the time 
of reception Ib- So, it will be very convenient to form the expression of the time transfer T{xa, xb) from a3^(t^) to 
XB{tB) in terms of the position of the receiver B at the time of emission Ia- For any quantity QB{t) defined along 
the world line of the station B, let us put Qb = Q{tA)- Thus we may write ^^(t^), ^^(t^), VB{tA), vb =| vb |, etc. 

Now, let us introduce the instantaneous coordinate distance Dab ~ xb — xa and its norm Dab- Since we want 
to know Ib — tA up to the order l/c**, we can use the Taylor expansion of Rab 



Rab = Dab + (^s - tA)vB + 7;{^b - tAf as + ^{tB - tA)^bB 



-{tB-tA)^aB + -{ 

where a^ is the acceleration of B and bs = das/dt. Using iteratively this expansion together with Eq. 



T{xA, Xb) == TixA, Xb) + —Dab ■ vb + it^Dab 



(Dab ■ Vb) 
^ab 



Dab ■ dB 



1 
ID 



(Dab ■ Vb) {v% + Dab ■ clb) + -^^ab ( 'Wb • as + ■:;Dab ■ bi 



ab 



^ vb['Tm(xa,xb) +Xj2{xa,xb)] 
c Dab 

— ^DabVb ■ [\m{xa,xb) + Xj2{xa,Xb)] 



we get 



(101) 



where T(xa, xb) is obtained by substituting xb, tb and Dab respectively for xb, tb and Rab into the time transfer 
function defined by Eqs. (88)-(92). This expression extends the previous formula [6] to the next order 1/c^. The 
second, the third and the fourth terms in Eq. (101) represent pure Sagnac terms of order 1/c^, 1/c"^ and l/c"*, 
respectively. The fifth and the sixth terms are contributions of the gravitational field mixed with the coordinate 
velocity of the receiving station. Since these last two terms are of order l/c**, they might be calculated for the 
arguments [xatXb)- 

VI. FREQUENCY SHIFT IN THE FIELD OF A ROTATING AXISYMMETRIC BODY 
A. General formulae up to the fourth order 



It is possible to derive the ratio qa/qb up to the order 1/c^ from our results in Sec. IV since Ia and Ib are given 
up to the order 1/c^ by Eqs. (93)-(96). Denoting by i'^")/c" the 0{n) terms in Z, qa/qb may be expanded as 



qA_ 
Qb 



= 1- 



1 Nab ■ [vA -Vb) , 1 



I -Nab 



Vb 



f(2) 7(2) 



VB 



7(3) 7(3) 

ly -VA -Ib -Vb 



+ -jNab ■ 
cr 



(4'^ • vb) [va - 2vb) + (J^A^ 



■VA Vb 



0(5). 



(102) 
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In order to be consistent with this expansion, we have to perform the calculation of u*^/u^ at the same level of 
approximation. For a clock delivering a proper time r, 1/vP is the ratio of the proper time cLt to the coordinate time 
dt. To reach the suitable accuracy, it is therefore necessary to take into account the terms of order \/c^ in goo- For 
the sake of simplicity, we shall henceforth confine ourselves to the fully conservative metric theories of gravity without 
preferred location effects, in which all the PPN parameters vanish except j3 and 7. Since the gravitational field is 
assumed to be stationary, the chosen coordinate system is then a standard post-Newtonian gauge and the metric 
reduces to its usual form 



Goo = l-^VK+^VK^ + 0(6), {GoO = ^^^^W + 0(5), Gy 



2--- W„,2 



1 + ^W ) 5^j + 0(4) 



where W given by Eq. (55) reduces to 

G [■ p*{x') 



W{x) = U{x) + ^ j ■ 



X ~ X' 



7+ 2 ' " 



(l-2/3)t/ + n + 374 



d^x' . 



(103) 



(104) 



and W is given by Eq. (74). As a consequence, for a clock moving with the coordinate velocity v, the quantity 1/vP 
is given by the formula 



I- — {W +-v^\ + 



from which it is easily deduced that 



P^-]W^ 



7 + i ) Wv"^ - \v^ + 2(7 + \)W ■ V 
'2 1 8 



+ 0(6), (105) 



^^l + M^--^- + r--r- 



1,2 1,2 



+ ^{(^ + l)(^-4i^^ - WbvI) + \{Wa- Wb) [Wa -Wb + 2(1 - I3)(Wa + Wb) +v\- 4] 



2{^ + 1){Wa-va-Wb-vb) 



'4 -^2 2 



0(6), 



It follows from Eq. (102) and Eq. (106) that the frequency shift Si^/i' is given by 



Siy va 



1 



V vb 
where {5v/v)(, is the special-relativistic Doppler effect 



(106) 



(107) 



(^ 



= --Nab ■ {vA -vb) + ^ 



2^A - 2^B " i^AB ■ {vA - vb)) (Nab ■ vb) 



{Nab ■ {vA - Vb)) ( x^a - ^"^ ^ {^^ab ■ vb)^ 



,4 1„,2„,2 



1 



g«A - I^aVb - ^Vb 

fl 1 

- {Nab ■ {vA - Vb)) {Nab ■ vb) ( ^v\ - -vl + {Nab ■ ^s)^ 



0(5) 



(108) 



and {5v)/v)g contains all the contribution of the gravitational field, eventually mixed with kinetic terms 



1 r. 



r(2) 



r(2) 



= ^{Wa~ Wb) -^ {Wa- Wb) {Nab ■ {va - vb)) - I'a' ■ va + I'b' ' vb 



+ ^ |(7 + ^){Wav\ - WbvI) + \{Wa- Wb) [Wa -Wb + 2(1 - I3){Wa + Wb)+v\-v 

-2 {Nab ■ {va - vb)) {Nab ■ vb)] + Nab ■ [(^b^ ■ vb) {va - 2vb) - 
l^^^ - 2(7 + 1)Wa) ■ VA - pj' - 2(7 + 1)Wb) ■vb} + 0(5) . 



(109) 
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It must be emphasized that the formulae Eqs. (105) and (106) are vahd within the PPN framework without adding 
special assumption, provided that (3 and 7 are the only non-vanishing post-Newtonian parameters. On the other 
hand, Eq. (109) is valid only for stationary gravitational fields. In the case of an axisymmetric rotating body, we 
shall obtain an approximate expression of the frequency shift by inserting the following developments in Eq. (109), 
yielded by Eqs. (95)-(100): 

Ia /c^ = -\m{xb,xa) - Xj^{xb,xa) H , ^^ /c^ = XsixB,XA) -\ , 

Z^Vc^ = \m{xa,xb) + X.hixA^xs) H , ^ /c"^ = \s{xa,xb) H , 

the function As being now given by Eq. (99) written with ai = 0. Let us recall that the symbol + • • • stands for the 
contributions of the higher multipole moments which are neglected. 

B. Application in the vicinity of the Earth 

In order to perform numerical estimates of the frequency shifts in the vicinity of the Earth, we suppose now that 
A is on board the International Space Station (ISS) orbiting at the altitude H = 400 km and that B is a terrestrial 
station. It will be the case for the ACES mission. We use vb = 6.37 x 10® m and m — r^ = 400 km. For the velocity 
of ISS, we take va = 7.7 x 10"^ m/s and for the terrestrial station, we have vb < 465 m/s. The other useful parameters 
concerning the Earth are: GM = 3.986 x 10^"* m^/s^, r^ = 6.378 x 10*^ m, J2 = 1.083 x 10"^; for n > 3, the muhipole 
moments J„ are in the order of 10^®. With these values, we get Wb/c^ ~ GM/c^vb — 6.95 x 10^^° and Wa/c^ ~ 
GM/c^TA = 6.54 X 10~^°. From these data, it is easy to deduce the following upper bounds: | Nab-va/c\< 2.6 x 10^^ 
for the satellite, | Nab ■ vb/c \< 1.6 x 10~® for the ground station and | Nab ■ {va — vb)/c |< 2.76 x 10^^ for the 
first-order Doppler term. 

Our purpose is to obtain correct estimates of the effects in Eq. (109) with are greater than or equal to 10^^^ for 
an axisymmetric model of the Earth. At this level of approximation, it is not sufficient to take into account the 
J2-terms in (Wa — Wb)/c^- First, the higher-multipole moments J3, J4, ... yield contribution of order 10~^^ in 
Wa/c^ ■ Second, owing to the irregularities in the distribution of masses, the expansion of the gcopotcntial in a series 
of spherical harmonics is probably not convergent at the surface of the Earth. For these reasons, we do not expand 
{Wa-Wb)/c^ inEq. (109). 

However, for the higher-order terms in Eq. (109), we can apply the explicit formulae obtained in the previous 
section. Indeed, since the difference between the geoid and the reference ellipsoid is less than 100 m, Wb/c^ may be 
written as [26] 

V.~ + ^^(l-3cos^^) + ^AI^., 
c-^ c'^rB Zc^r'^ & 

where the residual term AM^s/c^ is such that | AW^/c^ |< 10^^^. At a level of experimental uncertainty about 10~^*, 
this inequality allows to retain only the contributions due to M, Ji and S in the terms of orders 1/c'^ and l/c'*. As a 
consequence, the formula (109) reduces to 

^v\ 1 ,„. „, , 1 (5v\^'^^ 1 /(5i^^*^^ 



= (VK^ _ VKb) + - — 



/ g \/M \ / J2 

where the different terms involved in the r.h.s. are separately explicited and discussed in what follows 

)M 



By using Eq. (87), it is easy to see that {5v/v)\j is given by 



5v\^'^^ GM{rA 



-tb) 



M '^ATB 

Rab riA-VA+riB ■ vb 

+(7 + 1 



"^+1 '^A-rB\ . , ,, 

[Nab ■ (VA - vb)\ 



l + UA ■ TlB Ta + Tb , 

(111) 



ta+tb l + UA-riB 
The contribution of this term is bounded by 5 x 10^^** for 7 = 1, in accordance with a previous analysis [6] 



C. Influence of the quadrupole moment at the order 1/c 



(3) ■ 



It follows from Eqs. (98) and (109) that the term {Sv/vyj in Eq. (110) is given by 



(3) 

.12 



GM 



2rp 



J2 {Nab ■ (va-vb)) 



(7 + 1)GM 



(- + -) 

\rA TB } 



.h 



^j [3(fc.n,)^-l] 
1 



- W-r^B? 

TB 



1] 



{Nab ■ {vA - vb 
'■1 



TATB {l+UA- TIbY 

(k-riA + k -ubY frA_ rB_ 

1 + nA-n.B \rB ta 

rA{k-nBY +rB{k-nAf_\ fvA , tb 

TB 



1 



UA ■ riB 



TA +rB 



Rab , , (fc • n^ + fc • ub) 

-[UA -VA + riB ■ Vb^ 



TA 
2 



2 2 

l-UA-riB 



TA +rB 



TA TB 
TB TA 



2 TA 



-{nA-VA)[l-'Kk-nAf] 



2 rs 
Ral 



'-{nB-VB)[l-i{k-nBf] 



l + UA-riB 

rA+rB{2 + nA ■ ub) 

rA + rB 
rA{2 + nA -riB) +rB 



1 



-UA ■ riB 



rA +rB 



ua-va , Ub ■ Vb 



rA 



TB 



(k ■ nA){k ■ Ub) 



1 l-ik-UB? I, .l-{k-nA) 
-{uA- va) -(ns • Vb) 

2 rB 2 rA 



Ra 



B 



TA 

Rab 
Tb 



{k ■ Va) 

[k ■ Vb) 



, rA +rBi2 + nA ■ ub) , , 

k ■ ua h k ■ Ub 

ta + tb 

rA{2 + nA -tib) +rB 



k ■ UA + k ■ Ub 



rA + rB 
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(112) 



One has | va/c \~ 2.6 x 10^^, | vb/c \< 1.6 x 10^® and Kab = 3.77 x 10^^. A crude estimate can be obtained 
by neglecting in (112) the terms involving the scalar products ub ■ vb and k ■ vb- Since the orbit of ISS is almost 
circular, the scalar product ua ■ va can also be neglected. On these assumptions, we find for 7=1 



(3) 



< 1.3 X 10 



-16 



(113) 



As a consequence, it will probably be necessary to take into account the 0(3) contributions of J2 in the ACES 
mission. This conclusion is to be compared with the order of magnitude given in [6] without a detailed calculation. 
Of course, a better estimate might be found if the inclination i = 51.6 deg of the orbit with respect to the terrestrial 
equatorial plane and the latitude 7r/2 — 0b of the ground station were taken into account. 



D. Frequency shifts of order 1/c'' 



n(4) 



The term {5v /v) j^j in Eq. (110) is given by 



(4) 
M 



(7 + 1) 



V rA 
2 



GM 



rB 



GM{rA-rB), 2 2 



2rArB 



-K-4) 



+ \i^ [{rA-rB)' + 2{P-l){r\-rl)] 
2 \rArBj 



GM{rA+rB) 
rArB 

7+1 



2(7 + 1) _ rA -rB 
l + UA-riB rA + rB 
Rab 



[Nab ■ {va - vb)] {Nab ■ vb) 



l + UA ■ riB rA+rB 



{{uA ■ Va) {Nab ■ vb) - [Nab ■ {va - 2vB)]{nB ■ vb)} 



(114) 
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The dominant term {'y+l)GMv\/rA in Eq. (114) induces a correction to the frequency shift which amounts to 10""'^^. 
So, it wiU certainly be necessary to take this correction into account in experiments performed in the foreseeable future. 

The terms {Sh'/h')g is the contribution of the intrinsic angular momentum to the frequency shift. Substituting 
Eqs. (77) and (99) into Eq. (109), it may be seen that 



where 



(•^s)a = 



i^s)B 



X" J s 


= (Ts)a - i^s)B , 


(7 + l)'=f(l+''' 


i\ f kxUB TB 


3 ) ^ \l + nA-nB TA + TB 


k-{nAX ub) 
{l + UA- nsY 


'rA+TBi^ + nA-riB) 
ua + riB 

rA+TB 


Tb V r. 


3 \ f kxUA TA 


1 / ^ Xl + UA-riB TA+TB 


k ■ {uA X Ub) 

{l + UA -UbY 


TAi2 + nA -tib) +tb 
Ua + ■ riB 

TA +TB 



(115) 



k X UA 



(116) 



k X Ub 



(117) 



In order to make easier the discussion, it is useful to introduce the angle ip between xa and xb and the angle ip 
between the plane of the photon path and the equatorial plane. These angles are defined by 

cos ijj — Ua ■ riB , < -0 < TT , k ■ {ua x ub) ^ sin ip cos ip , < ip < tt . 

With these definitions, it is easily seen that 

k ■ (ua X Ub) V' 

= cos i„ tan — . 

I + ua-ub 2 

Let us apply our formulas to ISS. Due to the inequality vb/va < 6 x 10^^, the term {Ts)b in Eq. (115) may be 
neglected. From Eq. (116), it is easily deduced that 



|(^s)^|< (7 + 1) 



GS 



ta\ 2 + 3 |tanV'/2| 



TB 



1 1 + cos ij) I 



VA^ 



Assuming < ip < 7r/2, we have (2 + 3 | tan -0/2 |)/ 1 1 + cos-;/' |< 5. Inserting this inequality in the previous one and 
taking for the Earth GS/c^t\ == 3.15 x IQ-^^, we find 



1 f 5v 

A 



(4) 



< (7+1) X 10 



-19 



(118) 



Thus, we get an upper bound which is slightly greater than the one estimated by retaining only the term hoiV^ /c 
in Eq. (106). However, our formula confirms that the intrinsic angular momentum of the Earth will not affect the 
ACES experiment. 



VII. CONCLUSION 



It is clear that the world function fl(xA,XB) constitutes a powerful tool for determining the time delay and the 
frequency shift of electromagnetic signals in a weak gravitational field. The analytical derivations given here are 
obtained within the Nordtvedt-Will PPN formalism. We have found the general expression of n{xA,XB) up to the 
order l/c^. This result yields the expression of the time transfer functions Te{tA,XA,XB) and Tr(tB,XA,XB) up to 
the order 1/c^. We point out that 7 and ai are the only post-Newtonian parameters involved in the expressions of 
the world-function and of the time transfer functions within the limit of the considered approximation. 

We have treated in detail the case of an isolated, axisymmetric rotating body, assuming that the gravitational field 
is stationary and that the body is moving with a constant velocity Vr relative to the universe rest frame. We have 
given a systematic procedure for calculating the terms due to the multipole moments in the world function ^{xa, xb) 
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and in the single time transfer function T{xai xb)- These terms are obtained by straightforward differentiations of a 
kernel function. We have explicitly derived the contributions due to the mass M, to the quadrupole moment J2 and 
to the intrinsic angular momentum S of the rotating body. 

Assuming for the sake of simplicity that only [3 and 7 arc different from zero, we have determined the general 
expression of the frequency shift up to the order 1/c^. We have derived an explicit formula for the contributions of J2 
at the order \/c^ . Our method would give as well the quadrupole contribution at the order 1/c* in case of necessity. 
Furthermore, we have obtained a thorough expression for the contribution of the mass monopole at the fourth order, 
as well as the contribution of the intrinsic angular momentum S , which is also of order l/c"*. It must be pointed 
out that our calculations give also the vectors tangent to the light ray at the emission and reception points. So, our 
results could be used for determining the contributions of J2 and S to the deflection of light. 

On the assumption that the gravitational field is stationary, our formulae yield all the gravitational corrections to 
the frequency shifts up to 10^^* in the vicinity of the Earth. Numerically, the influence of the Earth quadrupole 
moment at the order l/c"^ is in the region of 10~^^ for a clock installed on board ISS and compared with a ground- 
clock. As a consequence, this effect will probably be observable during the ACES mission. We also note that the 
leading term in the fourth-order frequency shift due to the mass monopole is equal to 10~^^ for a clock installed on 
board ISS and compared with a ground-clock. As a consequence, this effect could be observable in the foreseeable 
future with atomic clocks using optical transitions. 
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